Direct Measurement of B(D° -> (pX°) and B(D + -> 
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The first measurement of B(D° —* 4>X°) and an upper 
limit for B(D + — > are determined from 22.3 pb^ 1 of 

e + e~ annihilation data at a C. M. energy of 4.03 GeV. The 
data were recorded by the Beijing Spectrometer (BES) at 
BEPC. A recoil charge method is applied for the first time 
to charm threshold data to determine the charge of the D 
meson in the recoil from 9054 ± 309 ± 416 reconstructed 
D°, D + mesons. The branching fractions B(D° -» (f>X°) = 
(1.71±S:n ± 0.17)%, and B{D+ -> 4>X+) < 1.8% are deter- 
mined from 10 events with a reconstructed D and a recoiling 
<f>. In addition, a 90% C.L. upper limit of B(D+ -» <j>e + X°) < 
1.6% is determined from a search for semileptonic decays of 
the D + . 



I. INTRODUCTION 

Much progress has been made in heavy flavor physics 
since the discoveries of open charm and bottom [EJ . How- 
ever, experimental results on inclusive decays of charm 
mesons are limited [^]||] . In an era of high precision ex- 
periments such as the B factories and the LHC, accurate 
measurements of 6-flavored particles can benefit from a 
better knowledge of charm decays and their branching 
fractions. In particular, the inclusive decay D — * <j)X [] 
has not been measured. This branching fraction serves as 
an independent check for the existence of additional ex- 



Throughout this paper, charge conjugation invariance is 
assumed, and charge conjugate states are included. 
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elusive decays of D mesons that contain a <f> meson. It is 
important for studies of gluonic penguin (b — > sg) decays 
of B mesons j|, and for time-dependent B®B S oscillation 
measurements Q that use a <j>£ pair to tag B® mesons. 
In addition, this branching fraction would be helpful in 
understanding the charm meson decay mechanisms. 

In this letter, we report a first measurement of the in- 
clusive <j) decay branching fractions of charged and neu- 
tral D mesons and a search for the exclusive semileptonic 
decay D + — > <f>e + X° based on an analysis of 22.3 p6 _1 
of data collected with the Beijing Spectrometer (BES) in 
e + e~ annihilations at y/s = 4.03 GeV. 



L>*(2007)°D*(2007)° 

and possibly DDn. The £>*(2010) _ can decay to either 
tt~ D or 7r°(7)-D - . The prompt or transition pion is not 
identified. Reconstructing a specific D meson does not 
necessarily determine whether the recoiling D meson is 
charged or neutral. In order to measure B(D° — ► 4>X°) 
and B(D + — ► <fiX + ) specifically, the numbers of neutral 
and charged D mesons recoiling against a reconstructed 
D meson, which tags the e + e _ — > DD events, and the 
types of the D meson from which the <fi mesons come, 
must be determined. 



II. THE BES DETECTOR 

The BES detector has been described in detail else- 
where 0. Here we briefly describe detector elements 
crucial to this measurement. 

BES is a conventional solenoidal detector operated at 
the Beijing Electron Positron Collider (BEPC) |§. A 
four-layer central drift chamber (CDC) surrounding the 
beam pipe is used for triggering purposes. A forty-layer 
main drift chamber (MDC) located just outside the CDC 
provides measurements of charged tracks and ionization 
energy loss (dE/dx ) with a solid angle coverage of 80% 
of 4-7T for charged tracks. A momentum resolution of 
1.7% \Jl + p 2 (p in GeV/c) and a dE/dx resolution of 
8.5% for Bhabha electrons are obtained for data taken at 
^fs = 4.03 GeV. An array of 48 scintillation counters sur- 
rounds the MDC and measures the time of flight (TOF) 
of charged tracks with a resolution of about 350 ps for 
Bhabha electrons and 450 ps for hadrons. Surrounding 
the TOF is a 12-radiation-length, lead-gas barrel shower 
counter (BSC) operated in limited streamer mode, which 
measures the energy of electrons and photons over 80% 
of 47r, with an energy resolution of a E /E = 0.22/ y/E(E 
in GeV), and spatial resolution of = 4.5 mrad and 
oz = 2 cm for electrons. Outside the BSC is a solenoidal 
magnet providing a 0.4 T magnetic field for the central 
tracking region of the detector. Three double-layers of 
muon counters instrument the magnet flux return, and 
serve to identify muons with transverse momenta greater 
than 500 MeV/c. They cover 68% of Att with a longitu- 
dinal (transverse) spatial resolution of 5 cm (3 cm) . 



III. THE RECOIL CHARGE METHOD FOR D 
TYPE IDENTIFICATION 



At yfs=4. 03 GeV charm mesons D° and D 
duced via the interactions 

e+e- -> D+D-,D°lf, 
D*{2010)+D- 
£>*(2010)+D*(2010) 



are pro- 
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FIG. 1. Monte Carlo Qrec distributions for (a) D + decays, 
and (b) D° decays. In both cases the D mesons are in the 
recoil against a fully reconstructed charm meson. 

A recoil charge method for identifying the type of the 
recoiling D meson has been devised for this measurement. 
At \fs =4.03 GeV, charged pions arising from Z?*(2010) + 
decay are very slow, and are mostly undetected in the 
BES detector. Only charged particles from decays of D + 
and D° are measured, and their total charge is correlated 
with the type of the mother D meson. Figures 1(a) and 
1 (b) show the Monte Carlo distributions of total detected 
recoil charge, Q re c, of D + and D° mesons against a fully 
reconstructed D meson, respectively. They are centered 
at +1 and 0, respectively, and have a spread of about 
+1 or -1. The recoil charge method selects neutral and 
charged D mesons according to 

Q rec = 0, or Q rec = F D = -1 for D° (1) 

and 



Qrec ■ F D < for 



(2) 



where Fd is the charm quantum number of the recon- 
structed D meson. The efficiency, e, and the misidentifi- 
cation rate, /, of the recoil charge method are obtained 
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from Monte Carlo simulations. For inclusive D decays, 
e and / are determined to be 0.74±0.02 and 0.25±0.02, 
respectively, and are the same for both charged and neu- 
tral recoiling D mesons. These values are confirmed 
with data using a kinematically selected data sample of 
e+e~ — > D + D~ and e+e~ — > D°D , where one D me- 
son has been reconstructed. The values from data, 
e=0.79±0.04, / = 0.16±0.04 agree well with the Monte 
Carlo estimates of e=0.795±0.015, / _= 0.204±0.015 
for the exclusive interaction e + e~ — > DD. The slightly 
higher efficiency for this sample is due to the absence of 
charged transition pions. 

When both a D meson and a recoil <j> are fully re- 
constructed, the efficiency of the recoil charge method 
is improved over that of the inclusive D events because 
there are fewer remaining tracks and they tend to clus- 
ter along the <f> direction and are therefore well within 
the detector acceptance. A Monte Carlo study of vari- 
ous D decay modes into final states containing a 4> has 
been performed, and the variations among their efficien- 
cies are included in the systematic errors. For events 
with a reconstructed 0, the recoil charge method identi- 
fies the recoil D meson type correctly (91±1±2)% of the 
time and misidcntifies a D in (9±1±2)% of the events. 
The first error is due to Monte Carlo statistics, and the 
second is systematic. 

To determine the numbers of neutral and charged D 
mesons on the recoil side, we use the relationships 
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(4) 



where No and N + are the events that pass the recoil 
charge tag as neutral and charged D meson candidates, 
respectively, N D o and N D + are the true numbers of neu- 
tral and charged D mesons in the recoil, and r + is the 
fraction of recoiling D + in the D tag sample. The values 
of N and N + are obtained directly from data. Using 
Eqs. (3) and (4) 7V D o and N D + in each D decay mode 
are determined. 



IV. DATA ANALYSIS 

A. Event selection 

Charged tracks are required to have good helix fits. 
These tracks must satisfy |cos6>| < 0.8, where is the 
polar angle, and must be consistent with originating from 
the primary event vertex. For charged particles, a com- 
bined confidence level calculated using the dE/drr and 
TOF measurements is required to be greater than 1% for 



the 7r hypothesis. For the kaon hypothesis, Lk > L n is 
required, where L is the likelihood for a particle type. 



B. Analysis of inclusive D meson events 

Charged and neutral D mesons are reconstructed in 
the £>° -» K-TT+, K-K-1T+TT+ and D+ -> K-n + ir+ de- 
cay modes. Figures 2(a), 2(b) and 2(c) show the mo- 
mentum distributions of selected K nir combinations with 
invariant masses within ±2.5 standard deviations of 
the nominal D° and D + masses. The lowest momen- 
tum peak in each figure corresponds to D mesons from 
e + e~ — > D* D , the middle peak is predominantly due to 
e + e~ — > DD , and the small enhancements at high mo- 
mentum are from direct e + e~ — ► DD production. 
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the selection criteria and the mass cuts mentioned in the text 



2. Momentum distribution of (a) the K it , (b) the 
A 7r + , and (c) the K~-k + -k + combinations that pass 



To reduce combinatorial backgrounds only D mesons 
in the lower two momentum regions are used in this mea- 
surement. Figures 3(a), 3(b) and 3(c) show the invari- 
ant mass distributions for selected Knir pairs. A binned 
maximum-likelihood fit to the distributions with a Gaus- 
sian signal function and a third order polynomial back- 
ground yields 10371 ± 357 events. Some D events en- 
ter the D signal region more than once due to K — it 
interchange and random combinations, and are doubly 
counted in the fit. The rates of double counting are eval- 
uated to be 14% for D° and 5% for D + , using data and 
Monte Carlo events. The largest variation in these rates 
between data and Monte Carlo simulation is found to be 
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3.6% and is included in the systematic error. Specific to 
the K~tt + tt + t:~ mode, partially reconstructed D° events 
combine with slow pions to produce an enhancement at 
the D° signal position with a resolution similar to that 
of the signal. A Monte Carlo simulation indicates this 
enhancement accounts for as much as 6.2% of the num- 
ber of D° — ► K ~ ir + ir + ir~~ events given by the fit. This 
variation is also included in the systematic error. 



An independent estimate of the number of recoiling D° 
and D + has been carried out using measured production 
cross sections of D*D and D*D and the branching frac- 
tions of D* and D mesons. The result is consistent with 
the recoil charge method. 

C. Measurement of branching fractions 

1. Reconstruction of <f> — > K + K~ events 
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FIG. 3. Invariant mass distribution of selected (a) K~tt + , 
(b) K~tt + it + -k~ , and (c) K~tt + -k + combinations. 

After subtracting the doubly counted events, a sample 
of 9054 ± 309 ± 416 tagged events is selected which con- 
tains in the recoil 6767±298±446 D° and 2287±83±155 
D + mesons, as determined by the recoil charge method 
using Eqs (3) and (4), where the first errors are statistical 
and the second systematic. The systematic errors are due 
to variations in the mass fits, the subtraction of doubly 
counted events, and the uncertainty of the efficiency of 
the recoil charge method. A summary of these events is 
presented in Table I. These recoiling D° and D + mesons 
are unbiased D decay samples for measuring branching 
fractions. 



TABLE I. Determination of D meson sample 



D tag 


number of 


r+ 


N D o 


N D+ 


type 


D events 








D" 


6895±288 


0.20 


5513 ±280 


1382 ± 69 


D+ 


2159±112 


0.42 


1254 ± 102 


905 ± 47 


Sum 


9054±309 




6767 ± 298 


2287 ± 83 



The 4> meson is reconstructed through its decay to 
K + K~ . Figure 4 shows the invariant mass distribution 
of K + K~ pairs. A fit of convoluted Breit-Wigner and 
Gaussian functions plus a third order polynomial back- 
ground gives a fit mass of 1.0194 ± 0.0002 GeV/c 2 and a 
total of 1108 ± 70 4> events. In this measurement, the (j> 
signal window is defined as the region from 1.00 to 1.04 
GeV/c 2 , as indicated by the arrows in Figure 4. 
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FIG. 4. Invariant mass of inclusive K + K~ pairs. 

2. Inclusive D — ► <pX 

Figures 5(a) and 5(b) show the mass of K + K~ pairs re- 
coiling against D° and D + candidates, respectively, and 
the full data are shown in Figure 5(c). 

To estimate the number of signal events, the K + K~ 
mass intervals 0.98-1.00 GeV/c 2 and 1.04-1.15 GeV/c 2 
are taken as background regions for the (f). The Knn 
mass regions from 1.7 to 2.1 GeV/c 2 , excluding regions 
within ±3ctm d of the fit D masses, are defined as side- 
band background control regions for the D mesons. As 
shown in Figure 5(c), 15 events are found as Dcj) candi- 
dates, and 14 events are selected as background sample 
outside the cf) mass region. Using the D side band events, 
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a total of 0.5±0.5 events has been estimated as the back- 
ground among the D candidates. Subtracting the back- 
ground contributions to both the D and the <fi, we obtain 
an excess of 10.2±4.0 events in the <f> signal region. 
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FIG. 5. Invariant mass distributions of K + K~ pairs recoil- 
ing against fully reconstructed (a) D°, (b) D + , and (c) D° or 
D+. 



events (Nd^)- 



B(D -> cAX) 



N 



D<j> 



N D x 



(5) 



where the <j> detection efficiency e^—0. 084±0. 006 as de- 
termined from Monte Carlo simulations. Assuming 
10.2±4.0 signal <j) events, the average branching fraction 
for the BES mixture of D° and D + mesons is measured 
to be 

B(D -» 4>X) = (1.34 ± 0.52 ± 0.12)%, 



where the first error is statistical and second systematic. 

Based on 9.7±lo D° -» 4>X° and 0.5±1;| D+ -» 
events, branching fraction and 90% C. L. limit 



B(D° -> <pT°) = (1.7118;?; ± 0.17)%, 



B(L>+ -» 0AT+) < 1.8%, 

are obtained, where the first error is statistical and the 
second systematic. The systematic errors arise from un- 
certainties in the numbers of D mesons and the error in 
the inclusive <j> efficiency. The final systematic errors are 
obtained by adding these uncertainties in quadrature. 
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FIG. 6. Invariant mass distributions of K + K~ pairs iden- 
tified as from (a) D°, and (b) D + . 

Figures 6(a) and 6(b) show the invariant mass of the 
K + K~ pairs from D + and D°, respectively, as identi- 
fied by the recoil charge criteria. After subtraction of 
backgrounds estimated using the 4> an d D side bands, 
9.7l£o D" -» <j)X° events and 0.5±^J D+ -> <j)X+ events 
are determined using Eqs. (3) and (4), and a multino- 
mial distribution that allows the fluctuations of the signal 
and the background. The branching fraction can be ob- 
tained from the number of D events (Nd) in the recoil 
against the tagged D mesons, and the number of D<j> 



3. Search for the decay D + — ► <f>e + X 

Of the 15 D(f> candidates selected four are accompanied 
by at least one charged track within | cos#| < 0.85. These 
tracks are tested against electron hypothesis where to 
identify electrons, a confidence level of greater than 1% 
is required, and L e > L v , using dE/dx information. None 
of the tracks are identified as electrons. 

TABLE II. Branching fractions and limits (90% C. L.) 



decay mode 


b. f. (%) 


Experiment 


D + -> <j>TT + 


0.61 ±0.06 


WA82 et. al. 


07T TV 


2.3 ± 1.0 


ACCMOR 


4>K + 


u.uio_ 019 


E687 


4>e + v 


< 2.09 


MK3 




< 3.72 


MK3 


CATT+Tr+TT- 


< 0.2 


E691 


4>e+X l> 
<j)X + 


< 1.6 

< 1.8 


this experiment 
this experiment 


D° - 


0.86 ±0.10 


ARGUS,CLEO, E687 




0.108 ±0.029 


ARGUS, CLEO, E687 




< 0.14 


ARGUS 


<t>n 


< 0.28 


ARGUS 


(j)LO 


< 0.21 


ARGUS 




1.7ll££ ±0.17 


this experiment 


D — > <f)X 


1.34 ±0.52 ±0.12 


this experiment 
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From no observed D + — > (pe + X event in a sample of 
2287 D + decays, and a detection efficiency of 0.0652 for 
the decay, a 90% C.L. upper limit of B(D+ -> <f>e+X°) < 
1.6% is placed. 

V. DISCUSSION 

These BES results, together with branching fractions 
listed in PDG98, are summarized in Table II. Compar- 
ing with existing measurements of exclusive D and D + 
decays containing a in the final states, as shown in 
Table II, these BES branching fraction values indicate 
little room for additional <j> decay modes of D° and D + 
mesons. 
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VI. CONCLUSION 



In summary, the absolute inclusive branching frac- 
tions of the D° and D + mesons decaying into a <f> 
have been directly measured. From a tagged sam- 
ple of 9054±309±416 DD pairs, 10.2±4.0 D -> <f>X 
events are observed, leading to the first measurement of 
B(D -> 4>X) = (1.34 ± 0.52 ± 0.12)% for a mixture of 
D° and D + mesons in the BES data sample, B(D° — > 
= (1.7ltjj-£ ± 0.17)%, B(D+ -» c^X+) < 1.8%, 
and B(D+ -» 4>e + X a ) < 1.6% at 90% C. L. 
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